arXiv:1503.01230v5 [hep-th] 16 Sep 2015 


UTHEP-669 


Matrix Geometry and Coherent States 


Goto Ishiku 

Graduate School of Pure and Applied Sciences, University of Tsukuba, 
Tsukuba, Ibaraki 305-8571, Japan 
Yukawa Institute for Theoretical Physics, Kyoto University 
Kyoto, 606-8502, Japan 


Abstract 

We propose a novel method of finding the classical limit of the matrix geometry. 
We define coherent states for a general matrix geometry described by a large-iV 
sequence of D Hermitian matrices Xf (y = 1,2,-•• ,D ) and construct a corre¬ 
sponding classical space as a set of all coherent states. When the classical space 
forms a smooth manifold, we also express various geometric objects on the classical 
space such as the metric, Levi-Civita connection, curvature and Poisson tensor, in 
terms of the matrix elements. This method provides a new class of observables in 
matrix models, which characterize geometric properties of matrix configurations. 
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1 Introduction 


Matrix models are conjectured to give nonperturbative formulations of string and M 
theories HIE]. In the matrix models, classical Riemannian geometry is replaced by a 
kind of quantum geometry described by matrices |3j and the matrix models are expected 
to realize a novel description of gravitational theories based on the matrix geometry. 
However, the relation between the matrix geometry and Riemannian geometry has not 
been fully understood so far @. Finding a clear relation between these two geometries, 
which look quite different at first sight, will help to understand the matrix models as 
theories of gravity. 

Although concrete examples of the matrix geometry (the matrix regularization) have 
been constructed for some manifolds !M> any algorithmic construction method has not 
been known for a general manifold. Also, its inverse problem, namely the problem of 
finding an associated classical (commutative) manifold for a given matrix geometry, has 
been only partly understood. For the latter problem, the Morse theoretic method [TO] 
and the method of taking an optimal gauge EHE2] have been proposed. These are useful 
to see some geometric aspects such as the topology or the shape of the classical space 
for a given matrix configuration. They are, however, not convenient enough to see more 
detailed information such as relations between matrix configurations and geometric ob¬ 
jects on the classical space such as the curvature tensors. For the case of two-dimensional 
surfaces, matrices which correspond to some geometric objects were explicitly constructed 
in [13] . Recently, another method of finding the classical space using a Dirac operator 
was proposed in [H10E]. This method is closely related to our proposal. 

In this paper, we propose a novel approach to the latter problem by generalizing the 
notion of coherent states. Coherent states are often considered for the Heisenberg algebra 
[p, q] — —ih ( see Appendix 0 and play an important role in taking the classical limit. 
They are approximate simultaneous eigenstates of p and q and can be defined as the states 
which saturate the uncertainty bound. They are considered as the quantum analogue of 
the points on the classical phase space. Indeed, there exists a one-to-one correspondence 
between the classical plane and a set of all coherent states. This means that, from the 
Heisenberg algebra, one can construct the classical geometry as a set of all coherent states. 

1 For infinitely large matrices, a very interesting solution to this problem was found in [4]. 
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The same construction is possible for the case of the fuzzy sphere, where the Bloch (spin) 
coherent states are in one-to-one correspondence with the points on S 2 (see Appendix iBl) . 
In this paper, we extend this construction to more general matrix geometries. 

For a Poisson manifold A4, its matrix regularization can be defined in terms of a large¬ 
ly sequence of matrices (for example, see [12]). Let Tjv be a linear map from C°°(Ai) 
to a set of all N x N matrices. When the sequence {T/vIat, where N is monotonically 
increasing to infinity, satisfies the four conditions shown in Eqs. (4.1)-(4.4) in [13] , {T/v}at 
is called the matrix regularization of Ai. The main property of the matrix regularization 
is the first two of the four, which take the following form: 

lim \\T N (f)T N (g) — T N (fg)\\ = 0, 

N—>oc 

Jim \\iN[T N (f),T N (g)]-T N ({f,g})\\ = 0. (1.1) 

N —>oo 

These conditions say that the map preserves the algebraic structure of functions and the 
Poisson bracket is well approximated by the commutators of matrices. In this paper, we 
consider the case where the manifold Ai is embedded in a flat space. For the embedding 
function : Ai —> R D , we denote by X M = T^{y^) the images of y A Then, the 
first condition in (11.11) means that the polynomials y^y v ■ ■ ■ are mapped to the matrix 
polynomials X^X V ■ ■ ■ up to the 1 /N corrections, while the second means that 

iN[X», X u ] - W^(X) 0, (1.2) 

where W^ w is the induced Poisson tensoj^]. Thus, the matrix regularization for an em¬ 
bedded space is given by a large-IV sequence of Hermitian matrices {jv which satisfies 
»• See [9] for more detailed description for the embedded case. In the matrix regular¬ 
ization of string and M theories, the matrices X' 1 are treated as the dynamical degrees of 
freedom and considered to realize a matrix regularization of the worldsheet of strings or 
the transverse worldvolume of membranes. 

In this paper, we will consider the inverse problem of the construction of the matrix 
regularization. Namely, we start from a large- N sequence of (bounded) matrices {X m }jv, 
which may or may not realize a matrix regularization of a certain manifold, and consider 

2 Of course, there is an ordering ambiguity in defining W with the matrix argument. However, ED 
also implies that X ^ commute with each other in the large-TV limit so that the ambiguity is negligibly 
small. 
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how a classical space M. can be associated with the matrix sequence. This kind of problem 
is relevant for some physically interesting situations, such as the emergent geometry in the 
gauge/gravity correspondence and the matrix realization of membranes. For this problem, 
(11.11) is not very useful since at the beginning we do not know what the space M. should 
be. Instead, we define coherent states for a given sequence {X^n as a generalization of 
the known cases and then define the classical space M. as a set of all the coherent states. 

For the cases where the classical space defined in this way forms a smooth manifold, 
we also express geometric objects on the classical space such as the metric, Levi-Civita 
connection, curvature and Poisson tensor, in terms of the matrix elements. If we consider 
a completely generic configuration of matrices, the classical space of our definition could 
be nonmanifold or an empty set. So in this sense, our result on the geometric objects 
cannot be applied to a generic matrix configuration whose classical geometries are not 
smooth manifolds. However, as will be discussed in the last section, our method is still 
valid for a small perturbation around matrices which have smooth classical geometry. 

The organization of this paper is as follows. In Sec. [21 we introduce the notion of 
coherent states for a general matrix geometry and define a classical space as a set of all 
coherent states. In Sec. El we express the metric, connection, curvature and Poisson tensor 
on the classical space in terms of matrix elements. In Sec. [H we show some examples of 
our construction. In Sec. 0 we consider the case of infinite-dimensional matrices. In Sec. 
El we summarize our result and discuss possible applications. 

2 Classical geometry of finite size matrices 

In this section, we consider the classical (commutative) limit of finite size matrices. We 
denote an index set by I := {a n |n G N}, where {a n } is a strictly monotonically increasing 
sequence of natural numbers. We assume that we are given a family of D matrices, 

,X^)\NeI}, (2.1) 

where xf/'- 1 (p = 1, 2, • • • , D) are N x N Hermitian matrices. We consider xj;, N> as 
regularized embedding functions with the target space R z5 , which possibly realize a matrix 
regularization of a certain manifold. We raise and lower the D-dimensional indices by 
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using the Kronecker delta (the flat metric on R D ), so that we do not distinguish upper 
and lower indices in the following. 

Let us define the position and its standard deviation in R /; for each statqj. Let PL be 
the N -dimensional Hilbert space on which the D matrices act and PL* be a subset of PL 
which consists of normalized state vectors: 


U* = {\a) G n\{a\a) = 1}. (2.2) 

For each |a) e PL*, we define the position of |a) in R 15 by 

x^(\a)) = (a\X^\a) (n = 1, 2, • • • , D) (2.3) 

and the standard deviations by 

^(H) = (<A x l\ a ) ~ {a\X„\a) 2 {n = 1, 2, • • • , D), 

D 

ff 2 (|a» = (2.4) 

fl=l 

Next, we introduce the notion of the coherent states. We recall that the canonical 
coherent states can be defined as the ground states of the Hamiltonians, H(p 0 , q 0 ) = 
~ Po ) 2 + |{q ~ qo) 2 , where po and go are real parameters (see Appendix [All . We 
generalize this definition. We first introduce the “Hamiltonian” 

= (2-5) 

1 

This is an N x N Hermitian matrix defined for each point y G R" 0 . The N x N identity 
matrix is omitted in the term of y fl . We denote the nth eigenstate and the eigenvalue of 
H{y) by | n,y) and Ejij): 

H(y)\n,y) = E n (y)\n,y) {n = 0,1, • • ■ , N - 1), (2.6) 

where we assume that the eigenvalues are ordered as E 0 (y) < Ei(y) < ■ ■ ■ < E N _ X (y) and 
the states are normalized as (n,y\m,y) = 5 mn . Our definition of coherent states is the 

3 In the following, we omit the superscript of (N) and make the IV-dependence implicit to avoid 
complexity. 
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following. We call |0,t/) a coherent state at y if and only if it satisfied 

lim E 0 (y) = 0. (2.7) 

N—too 

Onr definition of coherent states can be understood as follows. Suppose that the 
coherent state exists at a point y G R/°. Since the ground state energy can be written as 

EM = h 2 (|0.s)) + if>„(|0,s» - y„) 2 , (2.8) 

11=1 

Eq. (12.7P implies that cr(|0, y)) —> 0 and x M (|0, y)) —» y Thus, it follows that there exists 
a state whose wave packet is centered at y and shrinks to the point in the large -N limit. 
Conversely, suppose that there exists a state |a) G K* which satisfies cr(|a)) —> 0 and 
x^(\a)) —> yy. Then, ( a\H(y)\a ) goes to zero in the large-iV limit. Since E 0 (y ) can be 
written as 


E 0 (y) = min {ip\H(y)\ip), (2.9) 

\ip)£H* 

it also goes to zero in the large-iV limit; namely, a coherent state exists at y. Therefore, 
for each y G H D , if and only if a coherent state exists, there exists a wave packet which 
shrinks to a point in the large-iV limit. 

It is natural to define the classical space, which we denote by .M, as a subspace of H D 
on which there exist the shrinking wave functions. In terms of the coherent states, this is 
equivalent to the subspace of H D on which there exist the coherent states. Let us define 
a function / : H D —> R + as 


f(y) = lim E 0 (y). (2.10) 

N=> oo 

Then, we can write Ai as 

M = {y G R d | f(y) = 0}. (2.11) 

Note that one can compute f(y) from the given matrices in principle. So the expression 
(12.lip provides a relation between the classical space and the matrix configurations. 

4 Note that strictly speaking the coherent state defined here is not a single state vector but a set of 
the ground states which satisfy (1271) . 
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Let us comment on some properties of the function f(y). Suppose that E 0 (y) is a 
smooth function at a point y G R D , we can expand E 0 (y + e) in a Taylor series as 

Eo(y + e) = Eo(y) + e p d p Eo(y) + —e p 'e u d p d u Eo(y) + -e p e u e p d p d u d p Eo(y) + • • • . (2.12) 

Since E 0 (y + e) is the lowest eigenvalue of 

H(y + e) = H(y) + e„(T - X») + (2.13) 

we can compute the coefficients in (f? I ‘111 based on the perturbation theory ajj 


d P E 0 (y) =y l M-x li {\0,y)), 


N -1 


d„d v E Q {y) = 5 pu -2 J]Re 


77 = 1 


(0,y|X At |n,y)(n,?/|X i y|0,|/) 

E n (y)~E 0 (y) 


d,A;d p E 0 (y) = 

777,^0 72^0 


(0, y|A^|m, 1 /) (m, s/|X„|n, y) {n,y\X p \0, y) 
(E m (y)-E 0 (y))(E n (y)-E 0 (y)) 

(0, y\X u \n, y)(n, y\X p \0, y) 


-(O.bI^IO,!/) J2 
0 


(En(y) - E 0 (y)y 


+ ••• 


(2.14) 


where • • • represents the sum over all permutations of the indices /i, u, p. We can also 
expand the function f(y) as 

f(y + e) = f(y) + e p d p f (y) + ^e p e u d p d u f(y) + ^e p e u e p d p d u d p f(y) + • • • , (2.15) 


where d p f(y), d p d u f(y ) and d p d v d p f (y) are given by the large-iV limits of (12.14j) . Thus, 
we can compute the derivatives of f(y) as the corrections in the perturbation theory. 
This property plays an important role, when we define geometric objects in terms of the 
matrix elements in the next section. Note also that since x M (|0 ,y)) goes to y p for y e Ml 
as discussed above, the single derivative of f(y) is vanishing on Ml: 


d p f(y) = o (/i = l,2 ,■■■,£>, ye Ml). (2.16) 

Hence, f(y) is a quadratic function on a neighborhood of Ml. 

5 Here, we assume that the Hamiltonian is nondegenerate, but we can also treat degenerate cases in a 
similar way. 
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The coherent state at y is an approximate simultaneous eigenstate of the D matrices, 
where the eigenvalues are given by the components of y. In fact, it follows from (12.71) and 
the Cauchy-Schwarz inequality that, for any vector \ij}) G H*, 

lim — y M |0, y) — 0 (/x = l,2,■■■,£), y e M). (2.17) 

N—>o o 

From (12 .1 7j) . one can also show the following equations for y e M and for any state vector 

|V>) e U*: 

lim (ip\A(X - y)\0,y) = 0, (2.18) 

TV—>■ oo 

lim ((^\A(X)\0,y) - ^\0,y)(0,y\A(X)\0,y)) = 0, (2.19) 

TV—kx> 

lim {i/,\[A(X),B(X)]\0,y) = 0, (2.20) 

TV—>-oo 

where H(A) and B(X') are arbitrary polynomials of X M with finite (77-independent) 
degrees and coefficients. Here, we have assumed that A" M are bounded in the large-lV 

limit. 


3 Geometric objects 


In the following, we assume that there exists a neighborhood of A4 on which E 0 (y) is 
smootlo for any sufficiently large N. In this case, the classical space .Ad is a smooth 
submanifold of RT. In this section, under this assumption, we define various geometric 
objects such as the metric, Levi-Civita connection, curvature, and Poisson tensor on A4, 
in terms of the D matrices In the following, we assume for simplicity that the 

Hamiltonian (j2.5jl is nondegenerate on the neighborhood of A4, but the generalization to 
degenerate cases is straightforward. 


3.1 Metric 


We consider a real symmetric D x D matrix 


TV-1 


9nu(y) = 2 l r im E Re 

N—t OO z -' 


n=l 


(0, y\Xfj,\n, y)(n, y\X u \0, y) 
E n (y ) - Eo(y) 


(3.1) 


6 In order for the curvature tensors to be defined, Eo(y) needs to be at least 3 times differentiable on 
the neighborhood. 
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Note that this can also be written as 


9n*(y) = - d/j,d v f(y). (3.2) 

In the following, we show that the matrix g^ u (y) is a metric on AT 

Let us consider the expansion (j2.15(1 in which |e| = \ZS flu e tl e l/ is much smaller than the 
typical scale of AT If both y and y + e are contained in A4, we can regard e as a tangent 
vector on A4, and we denote e = ey. In this case, since f(y) = f(y + ey) = 0, we have 

9i*v(y)e K = l e !ll 2 - ( 3 - 3 ) 

On the other hand, let us consider the case in which e = e± is a normal vector at y E AT 
In this case, as shown in Appendix O f(y + ej_) is given by the distance between y and 
y + ej_ as 

f(y + ^) = \ M 2 - (3-4) 

By comparing this to Eq. (12.15ft with e = e±, we obtain 

9^{y)^±^L = 0. (3.5) 

Equations (13.31) and (13. 5 j) show that g^{y) is a projection to the tangent space TM. y at 
y G AT Thus, it satisfies 


g^{y)g u P (y) = gpp(y) (y e M). (3.6) 

Here, we again emphasize that the index v is raised using the Kronecker delta, so that 
g u p = g up . Since this implies that gp V {y) is positive and nondegenerate on A4, it gives a 
metric on AT The line clement on M. can be written as 

ds 2 = gp U (y)dy^dy u . (3.7) 


It is interesting that the positivity of the metric comes from the negativity of the second- 
order correction to the ground state energy in the perturbation theory of quantum me¬ 
chanics. 

Note that, for any y G M, the dimension of A4 is given by the trace of gpviyXX 


D 


dim M = '^Tgppiy). 

p=i 


(3.8) 


7 If At consists of some disconnected components, the trace of gpv(y) gives the dimension of the 
component which contains y. 












This expression and (13.ip relate the dimension of At with the matrix configurations. 


3.2 Levi-Civita connection 

Let A p (y) and B p (y) be tangent vector fields which satisfy g p v {y)A u {y) = A fl (y) and 
g p u (y)B u (y) = B p (y) for y e At. We define the covariant derivative on At by 

{\7 b AY = B v [d v A p + K p A p ). (3.9) 

The connection Y p p is chosen such that the image of the covariant derivative (13.91) is again 
a tangent vector, namely, 

g\(y B AY = (V b A)A (3.10) 

By using (13. 2ft and 

BP (y) {d P g^(y) - g^{y)d p g v a (y) - g^(y)d p g u p {y)} = o (ye At), (3.11) 

which is obtained by differentiating (13. 6 j) . we find a simple solution to (13.101) as 

r Y = (d°d^f)(d a d„d p f). (3.12) 

We also hnd that {X B g) P u = 0. This means that (13.121) is the Levi-Civita connection 
associated with the metric g pu . The expression (13.121) together with (j2. 14[) relates the 
connection and the matrix elements of X A 

3.3 Curvature 

Let A, B , C be tangent vector fields on At. The curvature tensor is defined by 

R{A, B)C = [Vn, V B ]C - V [AB] C, (3.13) 

where [A, B] represents the Lie bracket of the vector fields, 

[A, BY = A v d v B p - B v d v A A (3.14) 

Note that this is also a tangent vector satisfying g pu [A,BY = [A, BY- In terms of com¬ 
ponents, we can write the curvature as 

(R(A, B)cy = A v B p C a R p avp 
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= A"B P C° {(d p d„d x f)(d p d a d x f) - (d p d p d x f)(d„d a d x f)} . (3.15) 

By contracting the indices of R p C vp{y) using g pi/ (y), we can also write the Ricci tensor 
or the Ricci scalar in terms of the derivatives of f(y). These expressions and (12.141) give 
relations between the curvature tensors and the matrices. 


3.4 Poisson tensor 


We consider a real antisymmetric D x D matrix on A4 defined by 

W^(y) = lim c(0,y\[X' l ,X v ]\0,y), (3-16) 

JV—>-oo 

where c is a pure imaginary TV-dependent normalization constant. The constant c is chosen 
so that c(0, y\ [X^, X v ] |0, y) becomes O(N 0 ) in the large- TV limit (see the next section for 
concrete examples). W tiV {y) is a tangent bivector on Ai and satishes the Jacobi identity 
as we will show below. This means that W pu (y) is a Poisson tensor on A4. 

We first show that W pv (y) is a tangent bivector, namely, it satishes g p , P (y)W pu (y ) = 
W pi, (y) for y G A4. In order to prove this, the following relation is useful: 

c(0, y\[X p , X v ]\0, y){ 0, y\X„\n, y) ~ c(E n (y) - E 0 (y)){ 0, y\X p \n, y ), (3.17) 


where ~ stands for an equality for the leading-order terms in the large -N limit. This can 
be shown as follows. By using (12.191) and (12.201) . we obtain 

c(0,y\[X p y X^}\0,y)(0,y\X u \7i,y) ^ ^(0,y\{X^ -y^[X p : X u -y u \}\n,y). (3.18) 

Then noticing that {X u — y u , [X p , X v — y v ]\ = 2[X P , H(y)\, we obtain (13.171) . Now, let 
us calculate 


9^ l p(yW pv {y) = lim 

N—>oo 


2 Re 



(0,y\X p \n,y)(n,y\X p \0,y) 

E n {y) - Eo(y) 


(0,y\[X p ,X^]\0,y) 


(3.19) 


By using (13.17jl and the completeness relation |u, y){n, y\ = 1, we obtain 

g%{y)W p ^{y) = lim [2Re {c(0, y\X p X u \0, y) — c(0, j/|A^|0, y) (0, y|JW|0, y)}]. (3.20) 

N^-oo 


Since we have assumed that c is pure imaginary and X p are Hermitian, the second term 
is zero, while the first term is equal to W pu (y). Thus, we have shown that W pM (y) is a 
tangent bivector on A4. The Jacobi identity of the Poisson bracket can be shown in a 
similar way. We leave the proof to Appendix [D] 
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4 Examples 

In this section, we consider some examples. 


4.1 Fuzzy sphere 


Our first example is the fuzzy sphere, which is given by 


Xu = -= 


1 


L „ (yu = 1,2,3), 


(4.1) 


where L M are the iV-dimensional representation matrices of the SO (3) generators satisfying 
= ie^pLp and E^ =1 = ^ jT 1 lyvxiv- The normalization in (14.ip is chosen in 
such a way that the radius of the sphere becomes one: X ll X ll = InxN- The corresponding 
Hamiltonian (12.5ft is given by 


H(y) 


1 + \y | 2 2 Lpy» 

2 VN 2 - 1 ' 


(4.2) 


Here, we notice that for any SO(3) rotation, Lp —> A p v L u , there exists a unitary trans¬ 
formation which reproduces this rotation: 


/ iV ,,,/,- 1 = A / L „. 


(4.3) 


See Appendix | 

y = (y\y 2 ,y 3 ) 


for this transformation. Hence it suffices to consider the case where 
(0,0, r) with r — \y\ > 0. Then, the Hamiltonian reduces to 

1 + M 2 2L 3 \y\ 


H{y) = 


y/N 2 - l' 


(4.4) 


The ground state of this Hamiltonian is given by the highest eigenstate of L 3 with the 
eigenvalue (N — l)/2. Thus, we hnd that the function f(y) is given by 


m = 5(1 - Wf- 


(4.5) 


Therefore, the classical space A4, which is given as a set of zeros of f(y), is indeed a 
sphere with the unit radius. Here, we remark that E 0 (y ) is a smooth function for any 
point y apart from the origin, so that the smoothness condition assumed in the previous 
section is satisfied. 
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We can also compute the metric either from (13.ip or (13.21) . The result is given by 


9^(y) 



(4.6) 


For y G M, this is indeed a projection to the tangent space and (I3.7[) gives the standard 
line element on S 2 . It is an easy exercise to compute the curvature tensors from the 
definition in Sec. 13.31 and check that the Ricci scalar is equal to 2 and agrees with the 
known value for the unit sphere. 

The Poisson tensor (I3.16P is given by 


W^(y) = ie pvp lim 


2c 


TV— >• OO yj ^ 


(0, y\X p \Q, y). 


(4.7) 


If we put c = -y y/N 2 — 1, we obtain 


W^(y) = e^y p: 


(4.8) 


where we have used (12.171) . This is the standard Poisson tensor on S 2 embedded in R 3 . 


4.2 Fuzzy torus 

The fuzzy torus is defined in terms of the algebra of two unitary matrices U and V, 

VU = e ie UV, (4.9) 

where 9 = 2tt/N. The matrices U and V can be represented by the clock and shift 
matrices as 


U'rnn ^ j 

Fran dm+lm (4-10) 

where m,n = 0,1, • • • , N — 1 and 5 mn is the cyclic Kronecker delta which satisfies 5nu = 
Son■ By introducing four Hermitian matrices X p (ji = 1, 2, 3,4) as the real and imaginary 
parts of U and V as 

U = X 1 + iX 2 , V = X 3 + iX 4 , (4.11) 

we can regard the fuzzy torus as the noncommutative Clifford torus embedded in R 1 . Note 
that A" 1 and X 2 commute with each other and satisfy (A 1 ) 2 + (A" 2 ) 2 = Ijvxjv because of 
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the unitarity condition of U, and so do A" 3 and X 4 because of the unitarity condition of 
V. The Hamiltonian (12. 5 p is given by 

H(y) = \{U - z)(U f -z) + l(V- w)(V' - w), (4.12) 

where z = y 1 + iy 2 and w = y 3 + iy 4 . As shown in Appendix [Ej in this case the function 
f(y) is given by 

m = \( l-N ) 2 + i(l-M) 2 . (4.13) 

Hence, the classical space A4 is indeed the Clifford torus given by \z\ = |w| = 1. 

We can also obtain the metric from (13. 2 p as 

9ab{y ) = p - J, 9aa{y ) = 0, g a f){y) = (Sa/3 ~ > ( 4 ‘ 14 ) 

where a,b = 1,2 and a, ft — 3,4. We can see that g liu is a projection operator for y E M 
and gives the standard metric of the Clifford torus. 

In order to compute the Poisson tensor, we rewrite the algebra (14.91) using as 


[A x ± iX 2 , X 3 ± iX 4 ] - id(X 1 ± iX 2 )(X 3 ± iX 4 ), 

[W ± iX 2 , X 3 T iX 4 ] ~ -i9{X l ± iA 2 )(A 3 =F iX 4 ), (4.15) 


where we have neglected higher-order terms in 6. By putting c — i/9 in (13.16p . we obtain 
the Poisson tensor on the Clifford torus: 


W^(y) 


/ 0 0 

0 0 

2 4 14 

-y y y y 
\ y 2 y 3 -y 1 !/ 3 


y 2 y 4 -y 2 y 3 ^ 

- y l y 4 y 1 ^ 3 
o o 

o o / 


(4.16) 


5 Infinite-dimensional matrices 


As far as the Hamiltonian (12.5p has a discrete spectrum, it will be possible to apply the 
same formulation to matrices with an infinite size. In this case, we assume that we are 
given a one-parameter family of D operators xjp acting on an infinite-dimensional Hilbert 
space T~Lh, where h is a non-negative parameter. We replace the large -N limit considered 
in the previous sections with the limit of h — >■ 0 and the finite sums such as Yln=o wifi 1 
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the infinite sums appropriately. In this section, by considering a generic example, 

we will demonstrate that our formulation also works for infinite-dimensional matrices. 

In this section, we work with the notation as follows. We decompose the indices /i, is, p, 
which run from 1 to D , to two indices as p —> {A, /}. The indices A, B, C, ■ ■ ■ run from 
1 to 4, while /, J, K, ■ ■ ■ run from 5 to D. We further decompose A —> {a, a}, so that 
a,b, ■ ■ ■ run from 1 to 2 and a, (3, - ■ ■ run from 3 to 4. 

We consider a four dimensional noncommutative plane with some fluctuations embed¬ 
ded in R z;i , which is given by 

X A = q A + X A (A = 1,2, 3,4), 

X 1 = X T (/ = 5,6, - - - ,D). (5.1) 

Here q A are the operators satisfying 

[< q a , q b } = ie ab 9, [ q a , f] = ie^d', [q a , q a } = 0, (5.2) 

where e ab and e a/3 are 2x2 antisymmetric matrices with e 12 = 1 and e 34 = 1. 6 and 9' are 
the noncommutative parameters which are assumed to vanish in the classical limit h —> 0. 
We assume that the fluctuations {X IJ } are of the form 

X“= 1I a. = 1,2,.. • , D), (5.3) 

where k ■ x = Y^a=i ^ax a is the four-dimensional inner product and the Weyl symbols 
4>^{x) are real functions on R 1 . Here, we emphasize that (j)^{x ) are just C-number valued 
functions and they depend only on the four-dimensional coordinates xa- We assume that 
|0 M (x)| are small, so that we can treat perturbatively. The configuration (15.ip is often 
considered in the context of the description of gauge theories on the noncommutative 
plane in terms of matrix models |l6j . 

For the configuration (15.ip . the Hamiltonian (12.5p is given by 

H(y) = H 0 (y) + H 1 (y) + H 2 (y), (5.4) 

where H 0 (y),Hi(y) and H 2 (y ) are defined by 

H 0 (y) = ~ {q A ~ Va) 2 + \y], H x (y) = ]-{q A - y A ,X A } - X r y r , H 2 (y) = ^X 2 . (5.5) 
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The subscripts of the Hamiltonians indicate the degrees of the perturbation. The zeroth- 
order Hamiltonian H 0 (y) is just the two-dimensional harmonic oscillator. We introduce 
eigenstates of H 0 (y) by 

H 0 (y)\m,n,y) = Ej® n \m,n,y), 

E m,n = 9 (™ + ^ + d ' ( n + 0 + \y 2 i- ( 5 - 6 ) 

See Appendix[F]for explicit forms of the wave functions. In the following, in order to avoid 
the degeneration of Ho(y), we assume that 0 and O' are written as 6 = ah and O' = bh, 
where a and b are rational and irrational constants, respectively. Under this assumption, 
we can naively apply the formulas in the previous sections. 

In Appendix [Fl we compute the ground state energy of the full Hamiltonian H(y) up 
to the second-order perturbation in X fJ and derive the following form of f(y): 

f(y) = {:/ - <p’(y) + <P A (y)dA<P , (y)} h L ,(y) {/ - + <P B (y)B B <l> J (y)} + o(<#> 3 ), 

(5.7) 

where 

hij{y) = S u - d A <t> I {y)d A <t)j(y). (5.8) 

Hence, the classical space A4 is given by a four-dimensional surface in RU defined by 

yi = Mv) ~ <t> A (y)dA(t>i{y) + e>(</> 3 ) (j = 5,6, - - ■ ,d). (5.9) 

When = 0, the above equation reduces to y r — 0 and M. is just the four-dimensional 
flat space parametrized by {?m}- Equation (15.91) shows that Ad also fluctuates along the 
transverse directions when are turned on. 

From Eq. (I3.1j) or (13.21) . we can read off the metric on A4 as 

9ab{v) = Sab ~ 0A<My)<W 7 (y) + 0((p 3 ), 

9Ai(y ) = d A0i(y) - d A ((j) B (y)d B (j)i(y)) + C>(</> 3 ), 

9u{y) = d A (Pi(y)d A (pj(y ) + C>(0 3 ). (5.io) 

It is easy to check the projection property (j3.6[) for this metric. By substituting the metric 
(I5.10P into (13.71) . we obtain the line element 

ds 2 = (5 A b + d^ 1 (y)d B <j>i(y))dy A dy B + C>(</> 3 ), (5-11) 
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where we have used dy 1 = gi A (y)dy A obtained from (15.91) . By introducing new coordinates 
<j A by y A = a A + <p A (a), we can rewrite the line element as 

ds 2 = [(5% + d A (j) C {cr)){5 B c + d B <pc(^)) + d^^^dscp 1 (a)]da A dij B + O(cj) 3 ). (5.12) 

This is nothing but the induced metric associated with the embedding function given by 
X A (cr) = <j a + 4> a (ct) and X 1 (a) = 0 7 (a). 

Through a similar calculation, we can obtain the following form for the Poisson tensor, 

W AB (y) =B ab + B AC d c 0 B (y) - B BC d c <\>\y) + B CD d c 4> A (y)d D 4> B (y) 

- B AD (\F\y)d c d D <\) B [y ) + B bd ( f (y)d c d D (p A (y) + C>(0 3 ), 

W AI (y ) =B AB d B <t> I {y ) + B^ds^dc^iy) ~ B AD {y)d c d D cj ) I (y) + 0(0 3 ), 
W IJ (y) =B AB d A <j ) I (y)d B <t ) J (y) + O (0 3 ). (5.13) 

Here, we have defined a real antisymmetric matrix B ab by [ q A , q B ] = ihB AB and we have 
chosen the constant c in (j3. 16|) as c = —i/h. From (15.101) and (15.131) . one can easily check 
that W^ w (y) is a tangent bivector and its Poisson bracket satisfies the Jacobi identity. 

6 Summary and discussion 

In this paper, we proposed a novel method of finding the classical limit of matrix ge¬ 
ometries. We first introduced the notion of coherent states for general matrix geometry 
described by a large -N sequence of D Hermitian matrices (. X 1 ,X 2 ,--- ,X D ). We then 
defined the classical space AJ as a set of all coherent states. Assuming that the set of 
coherent states forms a smooth manifold and so does Ai, we also found expressions for 
various geometric objects on Ai (metric, connection, curvature and Poisson tensor) in 
terms of the D matrices. 

The usage of our result on the geometric objects is limited to the cases where the 
classical space is a smooth manifold. Probably, if we consider a general matrix sequence, 
the corresponding classical space is most likely to be an empty set. In order to have 
a nontrivial space, there must be a subspace of the Hilbert space on which the D ma¬ 
trices become mutually commuting with each other. Though this is not the case for a 
general sequence, there are physically relevant situations where such commuting matrices 
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naturally arise. In matrix models of the Yang-Mills type (i.e. models with commutator 
interactions), typical values of the commutators might become very small in the strong 
coupling region, since one can always take a normalization such that the coupling con¬ 
stant appears only in front of the commutators. As stated in S3. this gives a reasonable 
mechanism for the emergence of space in matrix models. Thus, in the strong coupling 
regime of matrix models, we can expect that nontrivial manifolds emerge and they can 
be visualized by using our method. 

Even if the space is nonempty, we needed to assume the smoothness condition to 
define the geometric objects. This condition is necessary to avoid some singular cases. 
The simplest singular example can be found for some special cases of diagonal matrices. 
Suppose that all of the D matrices are diagonal for any N. In this case, M. is just given 
by a discrete set of points corresponding to the positions of the eigenvalues in R D . Let 
us consider a further special case where the eigenvalue distribution becomes dense and 
continuous in a certain fixed region of R'° in the large -N limit, as the limit considered 
in the context of the large -N reduction [TBH2T] . In this limit, X4 becomes continuous 
manifold given by the eigenvalue distribution in the large -N limit. However, near the 
large -N limit, the ground state energy E 0 (y) is nondifferentiable almost everywhere on 
the dense region. Hence the formulation in Sec. [3] does not work. This example suggests 
that the off-diagonal elements are essential to describe the metric and the other geometric 
objects. Probably this kind of singularity should also be related with the topology change 
of membranes in M theory and we hope to find clear criteria to distinguish such singular 
configurations of matrices. 

Let us make some comments on our results. Firstly, our result can easily be applied to 
small perturbations around smooth configurations. For finite N matrices, let us consider 
an expansion of the matrices, Y M = X M + /r M ( X). Here, X^ are (a large-IV sequence of) 
matrices which satisfy the smoothness condition and have an associated classical manifold 
X4. hJ l (X ) are smooth polynomial functions which represent small fluctuations. In this 
case, by treating h^(X) as perturbation, from (]2.18j) - fl2.20|h one can show that the classical 
space for X^ is given by {y e R I) \y 1 '' = + (<5£ — v(y))h u (y ), y e M.}, where g /w is 
the metric of AT Thus, the perturbed geometry is again a smooth manifold. 

Secondly, if the matrices are thought of as the transverse coordinates of D-branes, the 
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classical space just corresponds to the classical shape of the D-branes. For example, the 
fuzzy sphere formed by matrices of DO-branes can be interpreted as the D2-brane arising 
via the Myers effect [22]. In this phenomena, the net D2-brane charge is zero, while there 
is a nontrivial gauge flux on D2-brane, which induces a coupling to the R-R 1-form. This 
gauge held will be realized as the Berry connection with respect to the Hamiltonian (12. 5 j) 
in our formulation. For example in the fuzzy sphere case, let us define A M i1|O), 
where | O) is the Bloch coherent state defined in (IB. 3D and n — 9, (j). Then we can find 
that the gauge held coincides with the Dirac monopole configuration, which agrees with 
the gauge held induced on D2-bra.nes [22]. Although we have discussed local properties 
of the coherent states in this paper, it will also be interesting to study global aspects in 
connection with the Berry phase of the Hamiltonian (12.51) . 

Now, let us discuss some possible applications of our result. Since all the geometric 
objects defined in this paper are invariant under the U(N ) gauge transformation X M —» 
UX^W , they can be regarded as observables in matrix models. These new observables 
will be useful in understanding geometric aspects of matrix models. For example, they 
can be used to hgure out the classical shape of D-branes in matrix models [12] and also 
to study the emergence of space-time in matrix models [23)21]. 

Our formulation also works for higher-dimensional gauge theories, and we can define 
the similar observables in gauge theories, which are useful in studying the gauge/gravity 
correspondence [25]. In the gauge/gravity correspondence, the background geometry of 
the string theory is expected to be emergent in the strong coupling limit of the corre¬ 
sponding gauge theory [26H28] . Though held configurations on the gauge theory side are 
considered to encode the information of the background geometry, a general prescription 
to see this is not known yet. Our observables will offer a new method of translating held 
configurations to the corresponding background geometry. 

Those observables will also be quite useful in computer simulations of matrix models 
[2T1I291I30] . By using a computer, one can easily compute the observables from given matrix 
configurations and see the geometric properties. Furthermore, although we considered the 
strict large- N limit when we defined the geometric objects, our formulation also works at 
finite N as an approximation, which is rather easily accessible from numerical simulation^. 

8 For example, the dimension of the fuzzy sphere (EH) with N = 50 can be read off from the value 
of (jfifi at the minimum of the ground state energy. One can find = 2.04, which gives a very good 
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Although we only considered the problem of the classical limit in this paper, fold¬ 
ing a general construction method of the matrix regularization is also a very important 
problem. In our formulation, coherent states introduce the notion of locality for matrix 
configurations. This idea may help to construct the matrix regularization for a general 
manifold. For instance, for a given manifold, one can first discretize the manifold to a foie 
lattice, and let the coordinates of each lattice cite to be diagonal elements of the matrices. 
In order to determine the off-diagonal elements, one may refer to the representation of 
the coordinate operators on the noncommutative plane in the basis of the coherent states. 
If any matrix geometry looks locally like the noncommutative plane, the off-diagonal ele¬ 
ments should be locally the same as those on the noncommutative plane. If so, by gluing 
those components along the tangent directions, it will be possible to construct the D 
Hermitian matrices which provide the matrix regularization. 

We pursue these directions and hope to report these issues in the near future. 
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A Canonical coherent states 

In this Appendix, we consider the canonical coherent states and summarize their proper¬ 
ties pTTj. We start with the algebra of the creation-annihilation operators, 


[a, a*] = 1. 


(A.l) 


The vacuum state is defined by 


a|0) = 0, 


(A.2) 


as well as (0|0) = 1. We consider the Hilbert space given by 


% = span{(a^) ?l |0 )\n G Z, n > 0}. 


(A.3) 


approximation of dimS 12 = 2. This kind of analysis will be possible for more general configurations. 
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The canonical coherent states are elements of T~L parametrized by a complex number 
z G C, and defined as 


| z) = e* at -*" a |0). 


(A.4) 


Let us introduce the Hamiltonian of the harmonic oscillator, 


H — h + -J , 

and its eigenstates 

H\n) = E n \n ), E n = K + 0 , | n) = -^=(a t ) n |0) (n = 1,2,- 

The coherent states can be written in terms of |n) as 

00 -,n 

\ z) = e -hw?y '\ n) _ 


The canonical coherent states have the following properties: 


Nil 2 l?2| 2 


(i) {Z]\z 2 ) = e “a "a + 5122 , 

f d 2 z 

(ii) / -|z)(z| = 1 (resolution of identity), 

J n 

(iii) a\z) = z\z), 

(iv) \z) is the ground state of the shifted Hamiltonian, 

H(z) = h - z)(a - z) + 0 , 

(v) | z) saturates the uncertainty inequality. 


(A.5) 


••)• (A.6) 


(A.7) 


(A.8) 

(A.9) 

(A.10) 

(A.ll) 

(A.12) 


The first and the second properties can be derived from (IA.7[) and the completeness 
relation Y2 n \ n )( n \ — 1- The integration measure in (IB.71) is defined as the flat measure 
for the real and imaginary parts of z, J d 2 z := fZm ez) d(Irri/~). The third, fourth 

and fifth properties follow from the fact that the unitary operator U(z) = e za] ~ za , which 
appears in Eq. (1A.4I) . is the translation operator W(z)aU(z) = a + z. 
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B Bloch coherent states 


Coherent states for the fuzzy sphere are called the Bloch (spin) coherent states |31‘ 9 
fuzzy sphere is described by the 577(2) Lie algebra: 


The 


[Z/j, Z/j] (B.l) 

We consider the spin J representation, so that the Li are (2 J + 1) x (2 J + 1) matrices. 
We introduce the standard basis for this representation space, which satisfy 


L 3 \ Jm) = m\Jm), 

L±\Jm) = \f(J =F m)(J ± m + 1)| Jnn ± 1), (B.2) 

where L± = Li±zL 2 , m = — J, —J+ 1, - • • , J and the states are normalized as ( Jm\Jnn') = 
5 m m'■ The Bloch coherent states are parametrized by the coordinates on S 2 , 72 = (0,0), 
where 0 < 0 < n and 0 < 0 < 2i r, and defined by 

|72) = e w* L —¥e-^ L +\ jj) (b.3) 


By using the Baker-Canbell-Hausdorff formula, they can also be written as 

172) = e * L - e - L 3i°s(i+\*\ a ) e -2 L +\ jj), 


where z = tan By acting the operators onto | JJ), one can rewrite them as 


172) = 


(1 + M 2 ) J 


E 

m=—J 


v J—m 


J 


E 


2 J 
J + m 


e 


2 J 

J + m 

J+m 


| Jm) 


cos 


6 


sm ■ 


J—m 




| Jm). 


m——J 

The Bloch coherent states have the following properties: 

8\ 09 J. \ 01 02 X 

_L nns — -I- p 1 ^ 2-^1) s in _ cin — 


(I) (72 1 |72 2 ) = ( cos — cos + e"^ ^ sm — sm 


(n) 


2J+ 1 

47T 


Is 2 


772172)(721 = 1 (resolution of identity), 


9 See also 


for the description of fuzzy sphere using coherent states. 


(B.4) 


(B.5) 


(B.6) 

(B.7) 
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(ID) y^L^Q) — J|?/||f2), where are given by (B.8) 

y 1 = \y\ sin0 cos 0, y 2 = \y\ sin0 sin 0, y 3 = \y\cosd, (B.9) 

(IV) 112) is the ground state of the Hamiltonian, 

' (R10) 
3 

(V) |f2) minimizes ^^(AL^) 2 , where Ais the standard deviation of L M . (B.ll) 

fi=i 


The first and second properties directly follow from (IB.51) . In property (H), the integra¬ 
tion measure is dehned as J di 2 = J^dO J^ w d(f> sin 0. The remaining properties follow 
from the fact that the unitary operator R(fi) := e h eeZrpL ~-h 6e zrpL + in Eq. (1B.3[) satishes 

R(n)L 3 R~\n) = 


C Derivation of equation ( 13.41) 


In this Appendix, we derive (13.4[) . As a preliminary step, we first consider the following 
object: 

{ 0 , yi \ H ( y )\ 0 , y 2 ) = ^0(2/2) + ^(2/2 - 2/ /i ) 2 ^ ( 0 ,2/i|0, y 2 ) + ( y 2 ~ 2/0(0, yi\X^ - y 2fl \0, y 2 ). 

(C.l) 

From (12 .171) . we obtain 


lim {y 2 -2/i)(0,J/i|0,J/2> = 0 ( 2 / 1 , 2/2 6 M). (C.2) 

TV—>-oo 

This implies that the coherent states at different points are orthogonal to each other: 

lim (0, 2 / 1 10, 2 / 2 ) = d [D) (yi - 2 / 2 ) ( 2 / 1 , 2/2 e A4). (C.3) 

TV->oo 

By applying (I2.17P and (IC.3P to the object (IC.ip . we obtain 

lim <0,2 /i|H(2/)|0,2/2> = \{y 2 - y^fd {D \yi - y 2 ) ( 2 / 1 , 2/2 e M). (C.4) 

iV—>-oo Z 

Now, let us compute 


f(y + e 1 ) = lim min 

V ' TV—>-00 |a)eW* 


(a|I/( 2 / + ex)|a), 


(C.5) 
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where y E M. and e± is a normal vector. Let us denote by 7-L co h the Hilbert space spanned 
by all the coherent states: 

'Hcoh = spanUO,?/)!?/ G M}. (C.6) 


The total Hilbert space can be written as 


= "Hcoh © a. 


(C.7) 


By definition, B is the Hilbert space which does not have any coherent state. We can 


show that when |e_i_| is sufficiently small, the minimum in flC.5|) is saturated by an element 
in 1-Lcoh- Hence, in the computation of (1C.51) . it suffices to consider the case where |ct) is 
an element of B co h- Any clement |cc) G "H C oh can be expanded by the coherent states as 


a) 


I 




(C.8) 


where a(y') is a function which vanishes unless y' G A4. Since |a) in (1C.51) is normalized 
as (ct| a) = 1, this gives a constraint on a(y'). By substituting the expansion (1C.81) into 
(1C.51) and using (1C.41) . we obtain 

f(y + e±) = ^ J d D y'\a(y')\ 2 (y^ - y 11 - e^) 2 . (C.9) 

Here, a(y') is the large -N limit of the function which saturates the minimum in (IC.5|) . 
It satisfies J d D y'\a(y')\ 2 = 1 and vanishes on the outside of M.. Obviously, the quantity 
on the right-hand side of flC.9jl is minimized when a(y') localizes at y. Thus, we finally 
obtain (j3.4[) . 


D Jacobi identity 

In this Appendix, we show that W^ u (y) defined in (j3.16[) satisfies 

W^{y)d, J W prT {y) + W^dyW^iv) + {y)d ll W vp {y) = 0. (D.l) 

This is equivalent to the Jacobi identity for the Poisson bracket defined by {A, B} = 
W^ u (d fl A)(d u B). We consider an arbitrary polynomial ^(A") and define a corresponding 
function by 

Hv) = j im (0,y|$(A)|0,?/). (D.2) 

N —^oo 
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We can show that 


lim c(0,y\[XrHX)]\0,y) = W^{y)d^{y). (D.3) 

TV—40o 


If we put $ = c[X tJ ',X u ] in the above equations, (ID.Hi immediately follows from the 
Jacobi identity of the matrix commutators. 

The relation (ID ,3|) can be shown as follows. The right-hand side of (ID. 3(1 is the large -N 
limit of 


c{0, y\{X“, JT]|0, y)d„(0, j/|4>(A')|0, y). 


(D.4) 


The derivative of |0,r/) can be read off from the formula in the perturbation theory as 

\n,y)(n,y\X u \0,y) 


&|0 ,y) = J2' 




E n (y) - E 0 (y) 


(D.5) 


By substituting this, and using (13 .1 7j) . we can see that the large-At limit of (1D.4|) is equal 
to the left-hand side of (ID.31) . 


E Classical geometry for fuzzy torus 


In this Appendix, we show that the function f(y) associated with the Hamiltonian (14.121) 
is given by (14.131) . First, note the following inequalities, 

Eo(y) = niin ty\H(y)\il)) 

> - min — z)(U * — z)\i[}) + - min — w)(V^ — w)\i[}) 

2 |' 2 \ip)eu* 




(E.l) 


The last inequality follows from the fact that the spectrum of any unitary matrix is 
contained in the unit circle on the complex plane. On the other hand, for any state vector 
| a) G "H*, we have 


E 0 (y) < (a\H(y)\a). 

Hence, if there exists a state vector |a) G H* which satisfies 

lim (a\H(y)\a) = ^(1 - \z\) 2 + ^(1 - |w|) 2 , 
iV— 4 oo Z 2 


(E.2) 


(E.3) 
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we can prove (14.13ft from the two inequalities (IE. 1)1 and (IE.21) . In the following, we 
explicitly construct a state |a) which satisfies (IE. 31) . 

We introduce the basis used in Eq. (14.101) . 


U\m) = e™\m), V\m) = | m + 1). 


(E.4) 


For convenience, we extend the range of m to a set of all integers by assuming the peri¬ 
odicity condition |m + N) = \m). With this notation, we introduce a state vector 


[f] 



e 


m=— 


JV-l ] 

2 J 


m 

N 


^S £ ) 2 e -*(arg w)m 


m), 


(E.5) 


where a is a constant of 0(V~N), c is a normalization constant determined by (a\a) = 1 
and [x] stands for the floor function defined by [x] = ma x{m G Z\m < x}. In the large -N 
limit, we can evaluate c as 


c 


2 



¥)‘~n dxe-<*-¥Y = 


where we have approximated the discrete sum by the integra 
calculations, we obtain 


0 


JV</f, ( E -6) 

V a 

By repeating similar 


lim (a\U\a) = e' arg2 , 

N^-oo 1 1 


lim (all/1a) 
N—>oo 


e ia,rgw 


This implies that |a) satisfies (IE.31) . 


(E.7) 


F Derivation of equation (15.71) 


In this Appendix, we derive (I5.7|) . Up to the second order of the perturbation, the ground 
state energy is given by 

Eo,o(y) = ~y~ + \ + °> y\ H ±(y) + H M 1 °, °> v) 

Z 1=5 

10 The integration range is naively given by [—1/2,1/2]. However, since we assumed that a = 0(y/N), 
the main contribution comes only from the neighborhood of the origin. Hence we can extend the range 
to [— 00 , 00 ] 
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E 

(m,n)^( 0,0) 


(0, 0, y\H 1 (y) I m, n, y)(m, n, y\Hi(y)\0, 0, y) 


£$(</) - ESM 


( 0 ) 


(F.l) 


We first construct the wave functions for the zeroth-order Hamiltonian Hq(ij) in 
Since H 0 (y ) is just the sum of two independent (shifted) harmonic oscillators, the wave 
function is given by the product of wave functions for each oscillator, 


^ m ,n(x,x;y) = 'SfV\x-,y\y 2 )ty { °'\x-,y 3 ,y A ). 


(F.2) 


'fm (i; y 1 , y 2 ) is the wave function of the mth excited state of the shifted harmonic oscil¬ 
lator, given by 

!/V) = c n h n (F.3) 

where c n = \j\J2 n 7i\\/Wrr and /r n (0 is the Hermite polynomial defined by h n (£) = 


^ 2 . The operators 

q A 

1 act on < F m , n (a; 

,x',y) 

(g 1 ^ 

,n) (%■> 

x; 

y) 

= X$ m , n (x, 

x;y), 


,n) j 

x; 

y) 

a 

(q 3 $m. 

,n ) 

x; 

y) 

= x$ m , n (x, 


(q 4 $m. 

,n ) j 

x; 

y) 

d 

= -M'-ir* 
ox 



(F.4) 


By using the wave functions (IF.21) . we can obtain 

(0, 0, y\e ik ' q \m, n, y) = ^^(iV9(h - tk 2 )) m (iV^(k 3 - ik 4 )) n e~^-^ +ik - y , 


(0, 0, y\{q a - y a , e lk ' q }\m, n, y) = i 
(0, 0, y\{q a - y a , e lk ' q }\m, n, y) = i 


k a 9 - 


kj - 


2mxa 
k\ — ik 2 
2 nxa 


(' 0,0,y\e lkq \m,n,y ), 

(I D,0,y\e ik - q \m,n,y ), (F.5) 


k 3 — ik 4 

where \a and Xa are constants given by xi — X3 — 1, X2 = Xa — —i- By using (IF.51) . 
we can compute the correction terms in (1F.1[) . In the classical limit where 9, 9' —> 0, the 
leading behaviors are given as follows: 

(0, 0, y\Hi(y)\0, 0, y) - y T (t>i(y), 

^=(0, 0,y\Hi(y)\l, 0, y) ~ -^= {My) ~ i<h(y) ~ V^diMv) ~ *Q 2 <Mj/))} , 
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=(o,o,j/|^i(j/)|o,i,j/) ~ -j= {Mv) -iMv) ~y\ d aMv) ~ id ^i{y ))}, 


(0,0,y\H 2 (y)\0,0,y) ~ ~(V{y)) 2 - 


(F.6) 


In deriving these, the following approximation formula for the delta function is useful: 


d k r -Tkl-^kl+ik-(y-x) _ \ r -Uya-Xa) 2 -jr(y a -x a ) 2 x(4)/ _ A 

(2vr) 4 ttW 1 V) - 


(F.7) 


By substituting (IF.61) into (IF.IF we finally obtain (15.71) . 
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